Kondo-resonance, Coulomb blockade, and Andreev transport through a quantum dot 
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We study resonant tunneling through an interacting quantum dot coupled to normal metallic and 
superconducting leads. We show that large Coulomb interaction gives rise to novel effects in Andreev 
transport. Adopting an exact relation for the Green's function, we find that at zero temperature, 
the linear response conductance is enhanced due to Kondo- Andreev resonance in the Kondo limit, 
while it is suppressed in the empty site limit. In the Coulomb blockaded region, on the other hand, 
the conductance is reduced more than the corresponding conductance with normal leads because 
large charging energy suppresses Andreev reflection. 

PACS numbers: 73.23.Hk, 73.23.-b, 74.80.Fp 
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Electronic transport of mesoscopic devices containing 
superconducting eloetrode has been an interesting sub- 
ject in recent yearsw Transmission of electrons through 
normal metal - superconductor (N-S) interfaces requires 
the conversion of normal current to supercurrent, which 
is called Andreev reflection.H With the recent advances of 
nanofabrication technique, quantum interference effects 
have been extensively studied in the-,mesoscopic N-S het- 
erostructures (see e.g. references inEl). In a phase coher- 
ent N-S structures, the phase of quasi-particles as well as 
Cooper pairs is preserved and transport properties de- 
pend strongly on the nature of quasiparticlqjphase. The- 
oretically, Landauer-Buttiker type formulaBu has been 
used extensively to describe quantum transport in many 
kind of N-S hybrid structures using non-interacting mod- 
els. (For review, see e.gJd'B.) 

Resonant tunneling through an interacting quantum 
dot (QD) or Anderson impurity has been intensively 
investigated recently. It has been shown that large 
Coulomb interaction gives rise to anomalous properties 
in transport. An example is the Kondo-resonant trans- 
port. Kondo-resonant transport has been, predicted the- 
oreticallyoia and verified experimentallyt3~E3 by conduc- 
tance measurements for artificially made Anderson impu- 
rities. On the other side, strong electron-electron interac- 
tions suppress conductance peaks where the systems are 
weakly coupled to the leads. It has been shown that 
the electron-electron interactions lead to conductance, 
suppression due to the orthogonality catastrophe.Ejt£l 
Stafford et alUB showed that the coherent transmission in 
artificial molecule structures is suppressed with increas- 
ing the system size by using the Hubbard-type model. 
Coulomb interaction has been found to P-lfY-a crucial role 
in the nature of the transmission phase ESiia It has been 
shown that Coulomb interactions give rise to anomalous 
effects in phase evolution through a quantum dot em- 
bedded in an arm of the Aharonov-Bohm interferometer, 
such as an inter-resonance phase drop. Non-equilibrium 
transport in an interacting quantum dot where both leads 
are superconductors has been studied recentlyJjpiising 
the non-equilibrium Green's function method Ej'E3 An- 
dreev reflection has been supposed to be negligible in the 



weak tunneling limit because large charging energy leads 
to Coulomb blockade of Andreev transport. In the mean- 
while, for a moderately coupled quantum dot, multiple 
Andreev reflections give rise to a novel subgap structure 
in the current-voltage curve due to resonant tunneling^, 
which are quite different from those of S-S contacts.Ej 
Resonant Andreev tunneling in strongly correlated quan- 
tum dot coupled to normal and superconducting leads 
has been investigated recently by Fazio and Raimondi.EJ 
Using the non-equilibrium Green's function formalism 
and equation of motion technique they have shown that 
the Kondo-resonant transmission is enhanced in the limit 
of large Coulomb repulsion due to the existence of a su- 
perconducting electrode. 

In this paper, we investigate coherent transport 
through an interacting quantum dot coupled to normal 
and superconducting electrodes based on the scatter- 
ing matrix formulation. We consider a model as shown 
schematically in Fig.l, where normal scattering and An- 
dreev reflection are decoupled. In the QD-N 2 boundary, 
only normal scattering is taken into account while An- 
dreev reflection is considered in the N2-S boundary. It 
is assumed that the N 2 -S boundary is perfect and the 
normal scattering doesn't occur at this boundary. This 
model is applicable to micro junctions where the length 
scale of nqwnal scattering and Andreev reflection is well 
separated.c3 With this model, Landauer type formula 
for the linear response conductance has been derived by 
BeenakkcrEa in the framework of non-interacting electron 
model. In the presence of interactions, this formula can- 
not be used in general because of the presence of inelas- 
tic processes. However, in the linear response regime 
(V = 0) with zero temperature, there is no phase space 
for inelastic processes and the formula can be equally ap- 
plied to the system containing interactions. The linear 
response conductance for the system under consideration 
can be written asE3: 
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where T n is transmission probability of nth channel. This 
equation is valid in the absence of an applied magnetic 



1 



field. We consider single channel case where the transmis- 
sion probability through the quantum dot is represented 
by T QD : 



Gns 



Ac: 
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The corresponding formula for the normal leads is the 
well known Landauer formulafl 
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Let's consider an Anderson impurity for the quantum 
dot with doubly degenerate level energy £o and on-site 
Coulomb repulsion strength U ~ e 2 /C, C being capaci- 
tance of the dot. At zero temperature the transmission 
probability Tqd can be obtained as follows iQwing to the 
fact that there are no inelastic processes.ESB Due to the 
absence of the inelastic scattering, the imaginary part 
of the self-energy for the Green's function at the Fermi 
energy ep is given by 



ImS(e F ) = -r/2, 



(4) 



where r = Tl + and T^/h and Tr/Ti are tunneling 
rate through left and right leads, respectively. With this 
condition the average occupation on the dot can be writ- 
ten as 



i) = -Im[logG r (e F )]. 



(5) 



At zero temperature, the transmission probability can be 
expressed in terms of the exact Green's function as 



T QD = T L T R \G r {e F tf 



(6) 



which leads to the final expression with the help of the 
Eq.(@) 



Tqd = p2 sin ip 



(7) 



where tp = n(n)/2. ((n) is average occupation of the 
dot.) 

Here (n) is calculated numerically by an equation of 
motion methodj 2 ^ which has been shown to be quite 
accurate for large U. We consider a symmetric cou- 
pling of the quantum dot to leads, that is Tl = Fr. 
From the calculated values of the average occupation, 
we display the conductances in the Fig. 2 as a function 
of Sf — Bq. The parameters used for calculations are 
U = 50r and W = 200r, with W being the band- 
width of the leads. Since the transmission probability 
reaches one for ef~Bq >> T, the conductance of normal- 
superconductor hybrid system goes to 4e 2 //i, which is 
twice of normal conductance. This is a result of per- 
fect transmission through the quantum dot in the Kondo 
limit. On the contrary, the conductances are suppressed 
in the empty site limit because of small transparency. 



Gns decays faster than Gn because transmission by An- 
dreev reflection requires two particle tunneling through 
the quantum dot. 

In real systems, nearly perfect Kondo-resonant trans- 
mission could not be realized though it is predicted by 
an exact relation at zero temperature. When the tem- 
perature is larger than the Kondo temperature, Kondo- 
resonant transmission does not occur and the transport is 
suppressed by Coulomb repulsion, rather than enhanced. 
In the case fcgT << T, thermal broadening can be ne- 
glected and the conductances can be obtained in the fol- 
lowing way with an approximate Green's function. If 
Kondo-like correlation is neglected, the transparency can 
be obtained by an approximate retarded Green's func- 
tioned which is similar to the Breit-Wigner type 
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Note that this equation coincides with the Breit-Wigner 
formula for U = 0. The self-consistent value of (n) is 
given by the relation 
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which leads to the expressiontS 
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where 



P 1 = iarctanfcM, 
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Then we can get the conductance through the Eq.(g) 

Fig. 3 displays the conductances obtained by the Eq.(||) 
and Eq.(j^). As one can see, Gns is suppressed more than 
Gn even in the "resonance" point. This phenomenon 
arises because Coulomb interactions in the dot suppress 
coherent transmission through the quantum dot. This 
would become a very general feature in transport through 
interacting system as far as the coupling to the leads 
are not so strong. If we consider larger system like a 
coupled chain of quantum dots, transmission probability 
will be more reduced due to orthogonality catastrophe. 
So one could say that in general Gns will be negligible 
compared to Gat in strongly interacting systems weakly 
coupled to leads. Normal conductance suppression with 
increasing the system size has been studied by Stafford 
et al.E-3 While the normal conductance is proportional to 
the transparency, Gns is second order of transmission 
probability. So Gns will decrease faster than Gn with 
increasing the system size. Even in the case of single 
quantum dot, we could see suppression of transmission 
due to the Coulomb interaction from our calculations. 
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For comparison, we plot the conductances of non- 
interacting case (U — 0) in the Fig. 4. As well known from 
the Breit-Wigner formula, one can see that Gns = 2Gjv 
in the resonance point because of perfect transmission. 
Comparing Fig. 3 and Fig. 4, one can conclude that the 
large charging energy suppress Andreev reflection even 
on resonance. 

In conclusion, we have discussed resonant tunneling 
through a strongly interacting quantum dot coupled to 
normal metallic and superconducting leads. We have 
found that in strongly interacting quantum dots, res- 
onant Andreev transport is qualitatively different from 
that of non-interacting system. Based on the scattering 
matrix formalism and adopting an exact relation for the 
Green's function, we have shown that at zero tempera- 
ture the linear response conductance is enhanced due to 
Kondo- Andreev resonance in the Kondo limit, while it 
is suppressed in the empty site limit. In the Coulomb 
blockaded region, on the other hand, the conductance is 
suppressed more than the corresponding normal conduc- 
tance even in the resonance point, because large charging 
energy suppresses Andreev reflection. 
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FIG. 1. Schematic diagram of the quantum dot(QD) coupled to normal(N) and superconducting(S) leads. In the N2-S 
interface, only Andreev reflection is considered. 



FIG. 2. Conductance Gns and Gjv obtained by the Eq.(7) and numerical calculation of cp for U = 50r and W = 200r. 



FIG. 3. Conductance Gns and G N obtained by the Eq.(8) with (6) for U = 50F. 



FIG. 4. Conductance Gns and Gn in the non-interacting (U = 0) limit. 
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